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We report ab initio calculations of the structural, electronic and magnetic properties of a graphene 
monolayer substitutionally doped with Co {Co sub) atoms. We focus in Co because among traditional 
ferromagnetic elements (Fe, Co and Ni), only Cosuh atoms induce spin-polarization in graphene. 
Our results show that, nearby the Fermi energy, the electronic structure of substitutional Co is 
equivalent to that of a vacancy in a tt- tight-binding model of graphene. Additionally, we investigate 
the magnetic coupUng between several Cosuh- We find that the total spin moment of arrays of Cosuh 
defects depends on the Co occupation of A and B graphene sublattices, in accord with the Lieb' s 
theorem for bipartite lattices. The magnetic couplings show more complex behaviors. 
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The peculiar electronic and magnetic properties of 
graphene monolayers have recently attracted much atten- 
New types of electronic and, particularly, spin- 



tion 



tronic devices based on graphene have been proposed. 
These graphene structures have also driven an increasing 
interest to study defects, which are always expected in 
real-life devices. Intrisic defects have been already widely 
studied^' ^i^iii^ and extrinsic defects, like substitutional 
atoms, are presently under intense research Re- 
cently, Gan et al}^ using high-resolution transmission 
electron microscopy (HRTEM) have observed substitu- 
tional Au and Pt atoms in graphene monolayers. They 
observed large activation energies for in-plane migration, 
which point towards strong covalent carbon-metal bond- 
ing and high stability for such substitutional defects. In 
spite of the growing interest about defects in graphene, 
little is known on magnetic properties, such as spin polar- 
izations and couplings, for the traditional bulk magnets 
Fe, Co and Ni placed substitutionally in graphene. 

In this work, we present a study of the structural, 
electronic and magnetic properties of substitutional Co 
{Cosuh) atoms in a graphene monolayer. We show that 
Co has spin polarization. This is in contrast to Fe and 
Ni substitutions which surprinsingly do not show spin 
polarization^^. We observe that there is a one-to-one 
correspondence between the expected behavior for single 
vacancies in a simple 7r-tight binding model of graphene 
and that found for the Coguh defects. The Co atom sta- 
bilizes a symmetric structure of the carbon vacancy. The 
electronic structure of the Cosuh defect at the Fermi en- 
ergy {E Fermi) IS dominated by a single level with a strong 
contribution from the Pz orbitals of the neighboring C 
atoms. Each Co^^^ defect shows a spin moment of 1.0 
jiB' The total spin, however, follows closely Lieb's theo- 
rem for bipartite lattices^^ and depends on the number of 
Co substitutions in each sublattice. Thus, magnetic cou- 
plings are predominantly (anti) ferromagnetic for Co sub 
defects sited in the (opposite) same sublattice. Our cal- 
culations also show the dependence of the couplings on 
the crystalline direction and relative position of the de- 



fects. 

We have used the SIESTA code^^iii for our calcula- 
tions. We have checked that the norm-conserving pseu- 
dopotential^- used for Co satisfactorily reproduces the 
band structure and spin moment of bulk Co. The pseu- 
dopotential includes nonlinear core corrections^^ for ex- 
change and correlation, with a pseudo-core matching ra- 
dius of 0.67 a.u. Computational parameters guarantee 
the convergence of relative energies for different magnetic 
configurations within the meV.^^ The results obtained 
with VASP code 20,21^ using a well converged plane- 
wave cutoff of 400 eV and the projected-augmented- 
wave scheme, are almost identical to those obtained with 
SIESTA. 

The relaxed geometry and the magnetization density 
for a Cosuh defect in a graphene monolayer is given in 
Fig. ma). The Co atom appears displaced 0.92 A from 
the carbon surface, and it has a symmetrical three-fold 
coordinated structure with Co-C distances of 1.77 A. The 
three C-Co-C bond angles are quite similar and have a 
value ~97^. The binding energy of the Co atom to the 
carbon vacancy in graphene is 7.6 eV. In comparison, the 
most stable adsorption configuration of Co on graphene 
has a much lower energy of 1.3 eV. Therefore, Cosuh de- 
fects are likely to be found in carbon nanostructures syn- 
thesized using Co-containing catalysts, as seems to be the 
case for Ni2^. Furthermore one could expect that Coguh 
defects can be fabricated in a controlled way by creat- 
ing vacancies in graphene and subsequently depositing 
Co. The presence of Coguh defects could be detected to- 
day by scanning tunneling microscopy (STM)^^, X-ray 
adsorption techniques'^ and HRTEM^^. 

Figure [1] (a) shows the spin-polarization pattern in- 
duced by the presence of a Cosuh impurity. The spin po- 
larization induced in the neighboring carbon atoms has 
a clear p^-like shape. The sign of the spin polarization 
follows the bipartite character of graphene: the polariza- 
tion aligns parallel (antiparallel) to the spin moment of 
the Co impurity for carbon atoms sitting in the oppo- 
site (same) sublattice than the Co atom. The total spin 




FIG. 1: (Color online) (a) Isosurface of the spin density in- 
duced by a Cosub defect. Positive and negative spin densi- 
ties correspond to light and dark surfaces with isovalues of 
zbO.008 e~/Bohr^, respectively. Panel (b) presents the spin- 
unpolarized band structure of an unreconstructed carbon va- 
cancy. Panel (c) and (d) show, respectively, the majority and 
minority spin band structure a Cosub defect in a similar cell. 
Filled symbols in panel (b) indicate the contribution of the 
Pz orbitals of the C atoms surrounding the vacancy, whereas 
empty symbols correspond to the sp^ character. In panels 
(c) and (d), filled and empty circles denote the contribution 
of hybridized Co 3d^2-C 2pz and Co 3d-C 2sp^ characters, 
respectively. Ef is set to zero. 



moment is 1.0 /j^b- Using Mulliken population analysis, 
0.44 jiiB are assigned to the Co atom; 0.20 /i^, to the 
three first carbon neighbors; and 0.36 jub, delocalized in 
the rest of the layer. The spin polarization decays with 
distance very slowly. 

To understand the origin of this spin polarization, we 
now analyze in detail the band structure of this system. 
Figures [1] (c) and (d) present the results for a Cosub de- 
fect in a 4x4 graphene supercell. Similar results are ob- 
tained using an 8x8 cell. For comparison, the panel (b) 
shows the spin-compensated band structure of a single 
unreconstructed carbon vacancy, the so called D^h va- 
cancy.^ For the Dsh vacancy, three defect states appear 
in a range of ~0.7 eV around the Fermi energy (E^?). 
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FIG. 2: (Color online) (a) Schematic representation of the 
geometry used to calculate the relative stability of ferromag- 
netic (FM), antiferromagnetic (AFM) and spin compensated 
(PM) solutions as a function of the positions of two Cosub 
impurities. Sublattices A and B are indicated by squares and 
triangles, respectively. One of the impurities is fixed in a cen- 
tral A-type site, whereas the other is moved along the (b) 
(n, n) and (c) (n, 0) directions. The empty circles represent 
spin compensated solutions and the full circles correspond to 
a fit with a Heisenberg model (see text). 



Two states appearing at 0.3 eV above at F have a 
large contribution from the sp^ lobes of the C atoms sur- 
rounding the vacancy. Other state at 0.35 eV below E^? 
at r shows a predominant Pz contribution. The Pz level 
corresponds with the defect state that appears pinned at 
Ei? for a vacancy in a simplified 7r-tight-binding descrip- 
tion of graphene. 

When a Co atom is bound to the vacancy, the defect 



states described above hybridize with the 3d states of Co. 
The two 2sp'^ defect bands, now an antibonding combi- 
nation of Co 3d and the original C 2sp^ vacancy levels, 
are pushed at higher energies, ~1.0 eV above Ef (see 
Fig. [11(c) and (d)). The singly occupied state, now hy- 
bridized mainly with the Co 3dz'2 orbital, remains pinned 
at the Ef and becomes almost fully spin-polarized. Thus, 
the Cosub impurities becomes analogous to a tt vacancy. 
The splitting between majority and minority Pz defect 
bands is ~ 0.50 eV. The spin splitting is much smaller 
for all other bands. 

We consider next the magnetic couplings between 
Cosub defects in a large 8x8 supercell with two Cosub 
impurities. We calculate the energy difference between 
different spin alignments as a function of the relative po- 
sition of the defects. Figure [2] shows the results along 
with a schematic representation of our notation. These 
energies are comparable to the spin-flip energies in fer- 
romagnetic transition metals Fe, Co and Ni. Several 
observations can be made: (i) when the impurities are 
located in the same sublattice (AA systems) the ferro- 
magnetic (FM) configuration is more stable than the an- 
tiferromagnetic (AFM) one; (n) if the Co atoms are in 
opposite sublattices (AB systems) it is very difficult to 
reach a FM solution, instead the system finds either 
a spin-compensated (PM) or an AFM solution; (iii) at 
short distances (< 3.0 A) the systems always converge to 
spin compensated solutions. 

In the FM cases, the total spin magnetization inte- 
grates to 2.00 /iB, with a spin population on every Co 
atom that remains almost constant ~0.50 /j^b (~0.30 /j^b 
for the three C-nearest neighbors). In other cases the 
total spin is zero. Thus, the total magnetic moment of 
the system follows the equation S = \^^b~ ^?ub\^ where 
^fub'' number of Cogub defects in the A(B) sub- 

lattices. Therefore, our total moment is consistent with 
Lieb' s theorem for bipartite lattices^^. This supports 
our analogy between the electronic structure of Co^^^ de- 
fects and single vacancies in a simplified 7r-tight-binding 
description of graphene. 

The spin magnetization density for some selected con- 
figurations is plotted in Fig. [3l Although the spin is 
fairly localized on the Co atom and the neighboring C 
atoms, the presence of the defect also causes a delocal- 
ized magnetization density with alternated signs in the 
two sublattices. The triangular pattern, that reflects the 
three-fold symmetry of the layer, shows different orienta- 
tions for A and B substitutions. This explains the very 
anisotropic AB interaction along the (n, n) direction seen 
in Fig.[2](b): the energy difference between AFM and PM 
solutions for (n, n)^^ configurations strongly depends on 
the relative position of the impurities, showing such di- 
rectionality. Similar patterns have already been observed 
experimentally^^'^^^^^'^^ in point defects (e.g. vacancies 
and hydrogens) on graphene by means of STM techniques 
and also theoretically discussed for tt- vacancies J^i^^j^ 
For Co sub in graphene, similar STM experiments should 
display the topology of the spin densities given in Fig. [3l 
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FIG. 3: (Color online) (a) Spin densities for configurations 
(a) (1, 1)^^, and (b) (-4, -4)^^ (see Fig.^a) for the nomen- 
clature). Positive and negative spin densities are indicated 
by light (gray) and dark (blue) isosurfaces corresponding to 
±0.001 e~/Bohr^, respectively. 

We have also investigated the magnetic interactions 
within the framework of a classical Heisenberg model, 
H = "^i^j JAA/AB{'^ij)^i^j^ where is the local mo- 
ment for a Co sub impurity at site i. The expression for 
the Tij dependence of the exchange has been taken from 
analytical RKKY coupling given in Ref. 31, except for 
the exponent of the distance decay, which is fitted to our 
ab initio results. The exchange interaction for AA sys- 
tems can be fitted with a l^ijl"^"^^ distance dependence 
(see the full circles in Fig|2] (b) and (c)). This distance 
dependence is in reasonable agreement with the \rij\~^ 
behavior obtained with analytical models for substitu- 
tional defects and voids^^^^^. In the case of AB systems 
a simple RKKY-like treatment fails to satisfactorily de- 
scribe the interactions, at least for the relatively short 
distances between defects considered in our calculations. 

Next, we explain the appearance of PM solutions in 
Fig. O The appreciable interaction between defect lev- 
els in neighboring impurities for AB systems opens a 
bonding- antibonding gap in the Pz defect band^^ and, 
thus, contributes to the stabilization of PM solutions. 
For A A systems, however, the bipartite character of the 
graphene lattice makes the interaction between defects 
much smaller. This explains the prevalence of solutions 
showing a local spin polarization for A A configurations. 
For very short distance between impurities, a larger 
defect-defect interaction opens a large gap and, in conse- 
quence, stabilizes PM configurations. It is interesting to 
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point out that similar behaviors were previously observed 
for vacancies in graphene and graphene ribbons described 
within a simplified 7r-tight-binding model. ^^'^^^^^ 

Similar behaviors to Cosub are expected for any other 
system that can be mapped into a simple tt- vacancy in 
graphene. This is the case of a hydrogen atom saturat- 
ing di C Pz state in graphene. ^-i^^ However, there are two 
important differences between these two systems. As Co 
prefers energetically to sit on carbon vacancies, Cosub 
defects could be fabricated by depositing Co atoms on 
areas where vacancies have been previously created by, 
for example, electron irradiation or Ar ion bombardment. 
The calculated diffusion barrier of hydrogen on graphene 
is 1.19 eV,^^ the experimentally observed activation en- 
ergy for the migration of substitutional metals, such as 
Au and Pt, in graphene is around 2.5 eV.^^. Thus, the 
diffusion is also probably much slower for Cosub defects 
than for hydrogen atoms. We expect that stable arrays 
of Cosub impurities could be fabricated and open a route 
to obtain in reality many of the interesting magnetic be- 



haviors obtained for 7r-models of defective graphene. 

In summary, only Co substitutional atoms in graphene 
respect to other traditional ferromagnetic materials, such 
as Fe and Ni, present spin polarization. The electronic 
structure of Cosub defects (see Fig. [1]) is dominated by an 
impurity level at with strong contribution from the 
Pz states of the neighboring C atoms, and it is a clear 
reminiscent of a single vacancy in a simplified 7r-tight- 
binding description of graphene. The study of the Co 
coupling in several sites demonstrates FM alignments for 
the AA cases, in contradistinction to AB cases which 
are either AFM or spin compensated. The AA cases 
were fitted to a Heisenberg model with a RKKY exchange 
constant. Our findings for Co^^^^ can be used to bring 
into contact the engineering of nanostructures with the 
results of 7r-models in defective graphene. 
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